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Dysfunction of the left ventricle (LV) weakens the cardiac function and 
affects the physical activity. Echocardiagraphy has been used to visualize 
the blood flow dynamics and to evaluate the cardiac function. 
However, the signal processing to suppress the clutter signals should be 
employed. In this study, we employed the singular value decomposition 
(SVD) clutter filtering to obtain the cardiac blood speckle images. We also 
employed the adaptive thresholding metric to determine the proper cutoff 
values at each phase during the cardiac cycle. Moreover, we employed 
a depth-dependent SVD clutter filter for more accurate estimation of 
the cardiac blood echo signals. The 2D blood flow velocity vectors were 
estimated by applying the block matching method to obtained blood speckle 
images. The obtained results show that the proposed filter suppressed 
the clutter signals from left ventricular wall significantly, 
and the contrast-to-noise ratio (CNR) was improved from -0.5 dB to 13.8 dB 
by the proposed SVD clutter filtering. 
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1. INTRODUCTION 

Echocardiography has been used to evaluate the cardiac function. The color-Doppler imaging 
is a credible technique for visualization of blood flow [1] and used for visualizing of the blood flowing in and 
out from the cardiac cavity. In echocardiography, the color-Doppler imaging is also employed to visualize 
an abnormal blood flow such as mitral regurgitation [2]. In the color-Doppler imaging, the flow velocity 
component only in the axial direction is estimated, and the velocity component parallel to the axial direction 
cannot be obtained. In this sense, the color-Doppler image doesn’t stand for the true velocities. 

The block matching method is able to estimate 2D velocities [3-6]. In this method, the motion of 
speckle patterns between successive frames is tracked and, thus, blood speckle images are required to obtain 
flow velocities. However, the intensity of the clutter echo is 40-60 dB higher than those originating from red 
blood cells [7]. Therefore, the signal processing to suppress the clutter signals should be employed. 
In our previous study, the high pass filter, which chooses the cut-off frequency adaptively based on 
the velocity of the heart wall [8], could visualize the blood B-mode images at a high-fra me rate [9, 10]. 

Recently, the singular value decomposition (SVD) gains attention in research fields of image 
processing as a low-lank approximation method. The SVD-based filtering was apphed to diagnostic imaging 
modahties [11, 12], and the SVD-based filter contributes to reduction of noise components contained in 
the diagnostic images. In the field of the medical ultrasound, the SVD-based method was employed as 
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a clutter filter [13-16] and its performance is better than that of the conventional finite impulse response 
(FIR) filter [13, 17]. The application of the SVD clutter filter is currently limited to blood flow in small 
vessels [13-15]. To assess the feasibility of the SVD filter in cardiac blood flow imaging, in this study, 
we applied the SVD clutter filter to echo signals the received from the human heart. In clutter filtering of 
ultrasound echoes from the heart waU, the heart wall exhibits remarkable movements during a cardiac cycle. 
For this reason, the fixed cutoff values cannot be adopted in the cardiac blood flow, and the cutoff values in 
a clutter filter should be altered during the cardiac cycle adaptively based on the information of the clutter 
motion as described in our previous study [8]. To solve such a problem, we employed the adaptive 
thresholding metrics in the SVD clutter filter. 

In the cardiac blood flow imaging, the optimum cutoff values were different at the positions of 
the echo data. Therefore, it is desirable to apply the clutter filters to spatially segmented ultrasound data. 
In this study, we employed a depth-dependent SVD clutter filter to obtain the cardiac blood echo signals. 
Also, an envelope-based block matching method was applied to the cardiac blood speckle images to obtain 
the 2D velocity vector maps. The performance of the proposed method was evaluated in the in vivo 
experiment on the human heart. 


2. RESEARCH METHOD 

2.1. Experimental setup 

A phased array probe at a center frequency of 3 MHz was used in the in vivo experiment. 
One diverging beam was transmitted using the transducer elements of the phased array probe [18, 19]. 
Element echo signals were received and applied the delay-and-sum (DAS) beamforming. The receiving focal 
points were aligned in the Cartesian coordinate to avoid a redundant effect due to non-constant lateral 
sampling interval in the polar coordinate [19]. As a result, beamformed RF signals along 241 scan lines were 
obtained. Under such conditions, a frame rate at 6250 Hz was achieved. 

2.2. SVD clutter filtering method 

A Let us define the beamformed RF signals s(x,z,t), where x,z, andt represent the lateral 
position, depth position, and time, respectively. The two-dimensional matrix S is created by rearranging 
the beamformed RF data s(.x,z,t) to the Casorati matrix with a size x 71^,71^), where 7i^,7i^,dind 7i^ 
is the number of pixels in the lateral and depth direction, and packet size of filter, respectively. By applying 
SVD, the matrix S can be decomposed into a production of three matrices as: 

5 = (1) 

where U and V are matrices composed of spatial and temporal singular vectors. 2 is a diagonal matrix 
composed of the singular values arranged in a descending order. The superscript represents transpose. 
The singular values obtained by such a procedure corresponds to energy of the received ultrasound signals, 
and the low-order, intermediate-order, and high-order singular components are assumed to be related to 
the clutter, blood, and noise signals, respectively. The clutter filter was accomplished by setting singular 
values corresponding to clutter signals to zero as: 

5 ' = ( 2 ) 

where I' is a diagonal matrix obtained by substituting singular values corresponding to clutter signal with 
zero. The filtered signals were obtained by reshaping S' to the original size. Figure 1 shows a schematic of 
the SVD clutter filtering. 

2.3. Thresholding of SVD clutter filter 

In this section, the thresholding metric for creating the matrix E' was introduced [15, 20]. 
In the SVD clutter filtering, the low-order components, which are assumed to be related to the clutter signals, 
were removed. Therefore, the boundaries between the clutter and blood signal components should 
be determined. As described in Section 1, in cardiac blood flow imaging, a cutoff threshold value of 
the clutter filter should be adaptively because the velocity of the heart wall changes rapidly during 
the cardiac cycle. 
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Eigure 1. Schematic of the SVD clutter filtering 


In this study, the blood subspace boundaries were determined based on the similarity of the spatial 
matrix U. The similarity between n-th and m-th column vectors of U was calculated as: 

C(n,m) ~ \U-n(x,z) I) x (|l/^(x,z) I - |[/^(x,z)|) (3) 

where U^(x,z) is the n-th spatial singular vector, and “ stands for averaging. The spatial singular vector 
corresponds to B-mode images. The two spatial singular bases, which belong to the same subspace, 
have high spatial correlation. Therefore, we determined the clutter and blood signal subspaces using 
the spatial similarity matrix C. Let us define the para meter matrix expressed as [15, 20]. 

fl if (n,m) G [l,a)^ 

^a,b = 11 if (n, m) G [a, b]^ (4) 

^ 0 otherwise 


The clutter and tissue and noise subspace boundaries, a and b, were searched so that the following 
normalized correlation coefficient Xn was maximized: 


Xn = 


x{c,aa b) 


Jx(C,d-z(o[a,b>^a,b) 


(5) 


where 


xiC,aa_b) = - \C(x,z)\) x (|a„j,Gc,z)| - |a„j,Gc,z)|) (6) 

Eigure 2 shows a schematic of an adaptive thresholding metric in this study. In the clutter subspace 
and in the blood subspace, the spatial bases in the same subspace were highly correlated. In Eigure 2, 
an area depicted as the red square was assigned to be the blood subspace. Then, the diagonal matrix I' was 
created so that the singular values except for the blood signal components were set at zero. 


Similarity of spatial vectors 
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Matrix for 2D fitting 




( 1 [1,a) 

1 \f in,rh^e [a,b] 
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Eigure 2. Schematic of the thresholding of SVD clutter filter in the case of a packet size of the filter of 500, 
variable a and b are chosen so that the similarity between C and ^ is maximized 
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2.4. Depth-dependent SVD clutter filter 

In Section 2.2, the SVD clutter filter is applied to the whole region of the received ultrasound data. 
However, in the cardiac blood flow imaging, the optimum cutoff values were different at the spatial positions 
of the echo data [8]. Therefore, it is desirable to apply the clutter filters to spatially segmented ultrasound 
data. In this study, we employed a depth-dependent SVD clutter filter to obtain the cardiac blood echo 
signals. As shown in Figure 3(b), regions of interest (ROI) were assigned by scanning a template window 
in the depth direction. The ultrasound data in the segmented regions were subjected to the SVD clutter 
filter and decomposed into a product of three matrix as: 

M; (7) 

where i = 1, 2 ^ the number of the ROIs. Blood signal components were extracted in each 

segmented data as: 

= (s) 

where M- is i-th the output signals of the SVD clutter filter. In this method, the blood signal components 
were estimated by the distribution of the singular values. To suppress the clutter signal components, 
we assigned the components with singular values of -55 dB lower than the maximum singular value as 
the blood signal components. Also, to suppress the noise components, the noise components were estimated 
using the slope of the distribution of the singular values [21]. 






Figure 3. Schematic of the SVD filter used in this study 


Finally, we concatenated output signals from each SVD filter to form one B-mode image. 
In this study, pixels in two adjacent ROIs were overlapped to reduce spatial discontinuity of the final output 
data. The numbers of overlapped pixels of adjacent ROIs were set to be 95% of a block size in the depth 
direction of ROI. In this case, one B-mode image was formed by concatenating 5% of output signals from 
each SVD filter. 

2.5. Estimation of the blood velocity vectors 

In this study, 2D velocity vectors were estimated by applying the block matching algorithm to 
the blood speckle images (envelope signals) obtained with the clutter filtering. The normalized correlation 
function between successive blood speckle frames was calculated with the kernel size of (6.0 mm x 7.9 mm) 
in the lateral and depth directions. The obtained normalized correlations were averaged temporally to 
suppress undesired components and interpolated with the reconstructive interpolation to obtain the subsample 
displacements [22-25]. Blood velocity vectors were obtained from the interpolated correlation function maps 
and displayed on the blood speckle B-mode images. 
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2.6. In vivo experiment 

The ultrasound echo signals from the left ventricular of a 27-year-old healthy volunteer were 
analyzed. The blood flow B -mode image were obtained by applying the SVD clutter filter to the beamformed 
RE signals. The B-mode image obtained with the SVD clutter filter was compared to that obtained with 
the EIR-based clutter filter, which chooses the cutoff frequency based on the velocity in the heart wall [8]. 
To evaluate the performance of the clutter filter, we measured the contrast-to-noise ratio (CNR) of 
the filtered images. CNR was expressed as: 

_ P^^Sblood)~P^i.Hissue^ 

std ^Stissue^ 

where s^iood •^tissue ^^ho signals which correspond to blood and clutter (background), and PWCs) 
denotes the power of the signal expressed as: 

PW(s) (10) 


Also, std(-) represents the standard deviation. 


3. RESULTS 

Eigure 4(a)-(c) show cardiac blood flow images at a systolic phase obtained without clutter filtering 
(original beamformed signals), with the EIR-based filtering, and with the SVD clutter filtering, respectively. 
In Eigure 4(c), the undesirable clutter signals at a depth of 100 mm were more suppressed than those shown 
in Eigure 4(b), whereas an appearance of the blood speckles were preserved. Regions of interest (ROIs) of 
the blood and clutter (background) were depicted in Eigure 4(a) as red rectangles. The CNRs of B-mode 
images in Eigure 4(b) and 4(c) became-5.8 and 12.3 dB. 



0 20 40 0 20 40 0 20 40 

Lateral position [mm] Lateral position [mm] Lateral position [mm] 


Eigure 4. B-mode image of left ventricular at a systole phase, (a) Original beamformed signals, 
(b) Cardiac blood flow images obtained with EIR-based filter, (c) With the SVD clutter filter, 
ROIs labeled with blood and tissue represent pixels corresponding to cardiac lumen 

and heart wall, respectively 


Eigure 5(a)-(c) show cardiac blood flow images at a mid-diastolic phase obtained without clutter 
filtering, with the EIR-based filtering, and with the SVD clutter filtering, respectively. In Eigure 5(c), 
signals which corresponding to the heart wall were also suppressed as compared with the B-mode image in 
Eigure 4(b). In this case, the blood speckle was less affected, but the SVD clutter filter was capable to 
suppress the clutter signals. The CNRs of B-mode images in Eigure 5(b) and 5(c) became -0.5 and 13.8 dB. 
Eigure 6(b) and 6(c) show cardiac blood flow images in systolic phase obtained by the SVD clutter filter 
without and with the segmentation of ROI in depth direction. 
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Figure 5. B-mode image of left ventricular at a mid-diastole phase, (a) Original beamformed signals, 
(b) Cardiac blood flows image obtained with FIR-based filter, (c) With the SVD clutter filter 



Lateral position [mmjLateral position [mm]Lateral position [mm] 

Figure 6. B-mode image of left ventricule in systole phase, (a) Original beamformed signals and obtained by 
SVD-based filter, (b) Without, (c) with the segmentation of ROI in depth direction 


Figure 7(b) and 7(c) show cardiac blood flow images in diastolic phase obtained by the SVD clutter 
filter without and with the segmentation of ROI in depth direction. In those results, is set to be 512. 
As shown in Figure 7(c), the contrast of echoes from blood cells is improved than Figure 7(b) by applying 
the depth-dependent SVD clutter filters to each ROIs. The CNRs of B-mode images in Figure 6(b) and 6(c) 
became 17.2 and 20.0 dB, and Figure 7(b) and 7(c) became 9.7 and 14.9 dB. Figure 8(a) and 8(b) show the 
2D blood velocity vectors of the cardiac blood flow at a mid-diastolic phase obtained with the FIR-based 
filtering and with SVD clutter filtering. In Figure 8(b), vortex flows in the left ventricle were captured weUin 
the proposed method. 
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Eigure 7. B-mode image of left ventricule in early-diastole phase, (a) Original beamformed signals and 
obtained by SVD-based filter, (b) Without, (c) With the segmentation of ROI in depth direction 
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Eigure 8. 2D blood velocity vectors on the blood speckle image obtained with, (a) EIR-based filter and with, 
(b) The SVD clutter filtering. Velocity vectors were obtained with the block matching technique 


4. CONCLUSION 

Visualization of the blood flow vector is useful to evaluate behavior of the cardiac blood flow. 
In this study, we applied the SVD-based clutter filter to the ultrasound echo signals from the human heart to 
obtain the cardiac blood flow image. We also employed the adaptive thresholding metric to determine 
the proper cutoff values through the cardiac cycle. Moreover, we employed a depth-dependent SVD clutter 
filter for more accurate estimation of the cardiac blood echo signals. The experimental results show that 
the proposed filter suppressed the clutter signals from left ventricular wall, and the CNR evalua ted using 
the blood speckle image was improved from -0.5 dB to 13.8 dB by the SVD clutter filtering. 
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